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Abstract 

The dynamical structure factor S(Q,to) of the S = 1 antiferromagnetic 
Heisenberg chain with length 20 at zero temperature is calculated. The lowest 
energy states have the delta-function peak at the region n > \Q\ > 0.37T. At 
\Q\ < 0.37T the lowest energy states are the lower-edge of the continuum of 
the scattering state, the strength of which decreases for large systems. This 
gives a reasonable explanation for the experimental fact that no clear peak is 
observed at the region Q < 0.3n. This situation is more apparent for valence- 
bond solid state. On the contrary for 5 = 1/2 antiferromagnetic Heisenberg 
chain the lowest energy states are always the edge of the continuum. 
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I. INTRODUCTION 



Haldane M first argued that integer-spin S antiferromagnetic Heisenberg chains have a 
singlet ground state which has a gap to a triplet excited states. This has now been well con- 
firmed both experimentally f^,[| and numerically 0-^] for S = 1 system. Recently, Ma et al 
carried out detailed inelastic neutron scattering experiments on Ni(C 2 H 8 N 2 )2N02(C104) 
(NENP), which is one of the most promising candidates for a Haldane gap system. They 
measured S(Q,u>) from Q = n down to Q = 0.3n where the intensity became weak. At all 
Q, the peak width is the resolution limit. Their experimental results coincide very well with 
numerical results of static correlation functions and low-lying excitations. The low-lying 
state gives the delta-function type peak in the region 7r > \Q\ > 0.3ir. It was stated that 
the gap near Q = is twice of that near Q = n [[| . This means that the lowest energy state 
near Q = is made up of two excitations near Q = n. Thus, near Q = the lowest energy 
state should be the lower edge of the continuum made by the scattering states. On the con- 
trary, however, near Q = it the lowest energy state has a strong delta-function peak. There 
must ,then, be a change in the property of the peak at a certain momentum. The Monte 
Carlo method is not appropriate to use to investigate the details of S(Q,u), so we used the 
diagonalization method and attempted to treat longer systems. This allows calculation not 
only of the first excited state but also of higher excited states and their contribution to the 
dynamical structure factor S(Q,u) [11]. The Hamiltonian we investigate in this paper is as 
follows: 

N 

H = Y: J[S?S? +1 + SfS? +1 + \S*S! +1 - /3(8t ■ S m ) 2 ] + D(S?) 2 , 

i=l 

Si = 2, Sat + i = Si. (1) 

The real NENP system is well described by A = 1,D = 0.18 J and (3 = 0. [0 In a 
previous paper || we calculated the cases A = 1,D = 0, ±0.2 J and /3 = 0uptoiV=18, 
and found that most of the spectral weight is concentrated in the lowest energy state. The 
ratio is more than 90% except near Q = 0. We then wanted to investigate what happens 
at small momenta and scattering intensity of higher excited states. The ground state was 



investigated by several authors fl3|J14}| about the point A = 1, D = 0, ft = — 1/3 . Following 



Affleck, Kennedy, Lieb and Tasaki [DJ we call this point the valence-bond solid (VBS) point, 
and we call the point X = l,D = 0,j3 = the antiferromagnetic Heisenberg (AFH) point. 
In this paper we carefully investigate the excitation spectra at these two points. 

We succeeded in diagonalizing N = 20 chain and determining the dynamical structure 
factor. The scattering intensity of the low-lying states decreases considerably in the region 
Q < 0.37r. Moreover, near 0.37T the second excited state approaches the first excited state. 
These facts mean that the single delta-function peak disappears at the region Q < Q c . 
This coincides with the experimental fact that Ma could not observe any clear peak in this 
momentum region. For comparison we calculate the dynamical structure factor at the VBS 
point ((3 = —1/3). Qualitatively the situation is very similar to the AFH case {(3 = 0). In 
VBS case the single delta-function peak is absorbed in the continuum at Q = OAtc. The 
lowest energy state becomes the lower edge of S(Q,u) and the intensity drops more rapidly 
and clearly. 

§2A summerizes our new method of representing states in the subspace of momentum Q. 
In our representation the Hamiltonian is always a real and symmetric matrix. State vectors 
have real elements. In the conventional representation Hamiltonian is hermitian and state 
vectors have complex elements. The memory necessary in our method is about half that of 
the conventional method. We can diagonalize N = 20 systems. 

§2B summarizes the method of Gagliano and Balseiro which is useful to calculate the 
dynamical structure factor. 

§3 compares the S = 1 antiferromagnetic Heisenberg (AFH) chain, S = 1 VBS chain 
and S = 1/2 AFH chain. We find that the S = 1 AFH chain is qualitatively very near the 
5 = 1 VBS chain, but the S — 1/2 AFH chain is completely different. We establish that 
the lowest energy state at Q > Q c has macroscopic strength and that at Q < Q c the lowest 
energy state is merely the lower edge of the continuum. The strength becomes microscopic. 
The lowest energy states of S = 1/2 AFH chain, on the other hand, are always the lower 
edge of the continuum and the strength is always microscopic. 



II. NUMERICAL METHOD 



A. Real symmetric matrix in momentum subspace 

As each site has three states, the Hamiltonian (1) is a matrix in 3 N dimensional space. 
Spin configuration of a state is represented by a trinary number with length N as 

| - 1,0,0,-1,1,0 >= |011021 3 > . (2) 

We introduce the shift operator T and inversion operator R as follows: 

T|011021 3 >= |110210 3 >, T N = I, 

R|011021 3 >= II2OIIO3 >, R 2 = /, RT' = T^R. (3) 

Both operators commute with the Hamiltonian (1). Here I is the identity operator. We 
can classify the subspace by S^ otal and momentum Q = 2n x integer / N . Each subspace is 
expressed by the following set of bases: 

N-l 

\a,Q >= const. ^ exp(-iQ/)T'|011021 3 >, T\a,Q >= exp(iQ)\a,Q > . (4) 
1=0 

We use the smallest trinary numbers in T'|trinary number> to represent a state \a, Q >. 
In this representation the Hamiltonian is a Hermite matrix which has the complex matrix 
elements, except the cases Q = 0, or n. The elements of eigenvectors are also complex. For 
N = 20, St otal = case, the dimensions of the subspace are 1.89 x 10 7 . If we can represent 
the Hamiltonian as a real symmetric matrix the memory and computing time should be 
reduced to half for the Hermite matrix. To find the real representation of Hamiltonian we 
classify the bases as follows. Let us call a base \a,Q > symmetric if 

R\a, Q >= exp(— iQl a )\a, — Q >, (5a) 

is satisfied for certain integer l a , and asymmetric if not. An asymmetric base \b,Q > must 
have its counterpart \b, Q > which satisfies: 

R\b, Q >= exp(-iQl b ) \b, -Q >, R\b, Q >= exp(-iQl b )\b, -Q >. (56) 
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For example, the state (011011) is symmetric and l a = 1. The state (011012) is asymmetric 
and (011021) is its counterpart and l b = 4. We generate a new set of bases as follows: 

| a', Q >= exp(iQl a /2)\a, Q >, 

\b',Q>= 2-^ 2 exp(tQl b /2)(\b,Q> +\b,Q >), 

\t, Q >= -12- 1 ' 2 exp(iQl b /2)(\b, Q>-\b,Q >). (6) 

In this new set of bases the Hamiltonian matrix is real and symmetric. Each column of 
the Hamiltonian has 2N non-zero off-diagonal elements at most. We therefore calculate 
the values and positions of non-zero off-diagonal elements and store them in an auxiliary 
memory before the Lanczos calculation. Thus only three vectors are needed in the core 
memory. 



B. Gagliano-Balseiro method 

We first calculate the lowest energy state |0 > in the subspace Sl otal = 0, Q = , namely 
the ground state using the Lanczos method. In the subspace S^ otal = 0,Q ^ we use 



Sq\0 > as the initial vector of the Lanczos calculation. Following Gagliano and Balseiro [IT 
we calculate the dynamical structure factors 

S\\(Q,u) = E I < n \ s Q\° > l 2 % - ( E n - Eo)/n\, 

n 

S x (Q,u) = J2\ <n\S Q \0> \ 2 5[u-(E n -E )/ti]. (7) 

n 

S\\(Q,uj) is given by the Green function: 

S\\(Q,w) = --ImG|,(Q,W), 

7T 

G,|(Q, z) =< 0\S z _ Q (z - H)~ l S z Q \U > . (8) 
For any operator A the Green function is given as a continued fraction by the Laczos method. 

„l Ai 1 A< n <0U t A|0> 

<0At ^A\0>= 1 1 — 2 , (9) 

z- ai 

Z — 0,2 _ _ 



\fo >= A\0 >, b = 0, 

|/ n+ i >= n\f n > -a n \f n > -6£|/ n _i >, 

a n =< fn\H\fn > I < fn\fn >, 
t>l+l =< fn+l\fn+l > I < fn\fn > ■ (10) 

We can determine poles and residues of this continued fraction and obtain the scattering 
intensities of each state. In the same way we can calculate S±(Q, u>). For high energy states 



this method does not give accurate results, but for low- lying states works very well [15 . 



III. RESULTS FOR SPIN CHAINS 

A. Isotropic Heisenberg chain 

In Fig.l we plot the scattering intensities of the lowest energy states of each momentum 
at isotropic Heisenberg point. At Q > 0.37T the ratio of scattering intensity of lowest energy 
state is more than 90% and not dependent on the size N. At Q < 0.3ir the ratio decreases 
considerably as N becomes large. The size dependence is strong and intensity seems to go 
to zero as N goes to infinity. In Fig. 2 we plot the first and second poles of momentum Q. 
We find that two states have almost the same energy at Q — 0.3tt, this is also evidence that 
the properties of the low-lying states change from this momentum. Figure 3 is the static 
structure factor S(Q). Figure 4 is the poles and residues of the dynamical structure factor 
S(Q,u). There is considerable distance between the first and the second excited state at 
Q > 0.3vr. 



B. VBS chain 

As the correlation length £ = 6.1 is comparable with size N = 20 the decrease in intensity 
of the lowest energy state is not clear for the AFH chain. But this situation becomes more 
clear if we calculate the ratio of the lowest energy state for the VBS chain. As shown in 
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Fig. 5 we find that the ratio goes to zero at Q < OAir. In the VBS case the correlation 
length £(= l/ln3) is much smaller than the system size and there is a clear decrease in 
the scattering intensity of the lowest energy state. In Fig.6 the first and the second excited 
energies are shown. Arovas, Auerbach and Haldane |nj obtained the structure factor S(Q) 
at the VBS point. Our result of S(Q) in Fig.7 reproduces their result. 

< |S: e S 5 |0 >= f- C ° S n - (11) 
1 Q Ql 5 + 3cosQ v ' 

Then, the variational energy becomes: 

1 <0|[^ Q ,[H,^]]|0> 10 

<^o = — 77— r- = — 5 + 3cosO . 12 

Q 2 <0\St Q S%\0> 27 v v; v ; 

The lowest energy state is found very near this variational energy at OAtt < Q < ir, but 
it enters into the continuum at Q < 0.47T. In this region S(Q,u) has only a vague peak 
near the variational energy, while the lowest energy state is far from this energy. There is 
considerable distance between the first and the second excited state at Q > OAn. 

C. S = 1/2 Heisenberg antiferromagnet 

It is known that the S = 1/2 Heisenberg model ( TC = JX^Sj • S i+ i, = 3/4) has the 
low-lying excitation e(Q) = |Jsin \ Q\. We calculate the ratio of the scattering intensity of 
this excitation in the states of momentum Q. As shown in Fig. 9. they decrease as N becomes 
large. Then, for this system the lowest energy state at fixed momentum is the bottom of 
the continuum made by scattering states of two excitations with spin 1/2 as was noted by 
Faddeev and Takhtajan [[Uj. This continuum has been justified experimentally [|17j]. The 



strength of each excitation is microscopic. The distance between the first and the second 
excitations decreases as the system enlarges. In Fig. 11 we show the structure factor of this 
system. It has a logarithmic peak at Q = it. In Fig. 12 we show the dynamical structure 
factor of N = 24 S = 1/2 AFH chain. 



7 



IV. CONCLUSION AND DISCUSSION 



From these numerical results we can conclude the following for S — 1 Haldane antiferro- 
magnets: 

1) At 7r > Q > Q c , the lowest energy state gives a scattering intensity of more than 90% of 
total intensity of the states of momentum Q. At Q < Q c the lowest energy state is the lower 
edge of the continuum and does not have the delta-function type peak. 

2) The value of Q c is about 0.37T for isotropic Heisenberg point and about OAn for VBS point. 

It is very interesting that the kink in e(Q) at Q = Q c and (3 = is too weak to be 
observed. 

The author is grateful to Ian Affleck and Stephen Haas for stimulating discussions. Most 
numerical calculations were done by HITAC S-3800 at the computer center of the University 
of Tokyo. This work was supported in part by Grants-in-Aid for Scientific Research on 
Priority Areas, " Computational Physics as a New Frontier in Condensed Matter Research" 
(Area No 217) and "Molecular Magnetism" (Area No 228), from the Ministry of Education, 
Science and Culture, Japan. 

APPENDIX: PROOF OF REAL-SYMMETRINESS OF HAMILTONIAN 

We show that Hamiltonian is real and symmetric in the bases given in equation (6). It is 
sufficient if we prove the symmetriness, because the Hamiltonian matrix is Hermitian. We 
first consider matrix elements between two symmetric states \a[,Q > and \a' 2 ,Q >'■ 



< a' 2) Q\n\a[, Q >= exp(iQ(h - Z 2 )/2) < a 2 , Q\H\a u Q > . 



(Al) 



Using Eq. (4a) and H = R7iR we have 



(Al) = exp(-iQ(Zi - Z 2 )/2) < a 2 , -Q\H\a^ -Q > . 



(A2) 



As the Hamiltonian is Hermitian in conventional representation we have: 



(A2) = expHQ^ - Z 2 )/2)< a ± , -Q\H\a 2 , -Q > 
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= exp(-iQ(h - l 2 )/2) < ai, Q\H\a 2 , Q >=< a[, Q\H\a' 2 , Q > . (A3) 

Thus we know that < a' 2 , Q\Ti.\a[, Q > and < a[, Q\H\a 2 , Q > is the same and is real. 
Next we consider elements between symmetric and asymmetric states: 

< b', Q\H\a, Q >= exp(iQ(l a - l b ) /2)2~ 1 / 2 {< b,Q\+ < b, Q\}H\a, Q > . {AA) 

Using Eq.(4b) and H = RHR we have 

{AA) = 2- 1 ' 2 exp(-iQ(l a - k)/2){< b, -Q\+ < b, -Q\}H\a, -Q > . (A5) 

By Hermiticity we have: 

(A5) = 2~ 1 ' 2 exp(-iQ(l a - l b )/2)< a, -Q\H{\b, -Q > +\b, -Q >} 

= 2" 1 / 2 exp(-^(/ a - l b )/2) < a, Q\H{\b, Q > +\b, Q >} =< a', Q\H\b', Q > . (A6) 
In the same way we can prove the symmetriness of elements between \b\ Q > and \a', Q >. 

<b',Q\H\a',Q>=iexp(iQ(l a -l b )/2)2' 1 / 2 {< b,Q\- <b,Q\}H\a,Q > 

= -12- 1 ' 2 exp(-iQ(l a - l b )/2){< b, -Q\-<b, -Q\}H\a,-Q > 

= -12- 1 ' 2 exp(-iQ(l a - h)/2)< a, -Q\H{\b, -Q > -\b, -Q >} 

= -i2' 1 ' 2 exp(-iQ{l a - l b )/2) < a, Q\H{\b, Q>-\b,Q>} =< a', Q\H\b\ Q > . (A7) 

In the last stage we must prove the symmetriness of elements between states b[,b 1 and 
states b 2 ,b' 2 . We have three cases which can be proved in the same way: 

<b' 1 ,Q\H\b' 2 ,Q>=< b' 2 ,Q\H\b' 1 ,Q>, 

<b[,Q\H\b' 2 ,Q>=< b' 2 ,Q\n\b[,Q>, 

< b[, Q\H\b' 2 , Q >=< 4 Q\n\b[, Q > . (A8) 
Thus, in our representation the Hamiltonian matrix is symmetric and real for any momentum 

Q. 
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FIGURES 

FIG. 1. Fraction of the total scattering intensity concentrated in the lowest energy state for 
N = 16, 18, 20 chains for isotropic Heisenberg point. In the region Q > 0.37T the ratio is more than 
90% and a stable function of momentum Q, but it decreases considerably as N increases in the 
region Q < 0.3-ir. 

FIG. 2. Energy of the first and second excitations of 5 = 1 isotropic Heisenberg point. They 
are very close together at Q = 0.3-7T. 

FIG. 3. Structure factor S(Q) of 5 = 1 isotropic Heisenberg model. It has a Lorentzian peak 
at Q = it and behaves as Q 2 near Q = 0. 

FIG. 4. Dynamical Structure factor S(Q,u) of N = 20, 5 = 1 isotropic Heisenberg model. 
It is made of delta-function peaks at each excited state. The strength is expressed as a percentage 
of total strength at the momentum Q. 

FIG. 5. Fraction of the total scattering intensity concentrated in the lowest energy state for 
N = 16, 18, 20 chains for VBS point. In the region Q > OAir the ratio is more than 90% and a 
stable function of momentum Q. In the region Q < OAir it decreases considerably as N increases. 

FIG. 6. Energy of the first and the second excitations of 5 = 1 VBS model. They are very 
close together at Q = OAir. 

FIG. 7. Structure factor S(Q) of 5 = 1 VBS model. It has a Lorentzian peak at Q = tt and 
behaves as Q 2 near Q = 0. It obeys eqn.(ll) exactly. 

FIG. 8. Dynamical structure factor S(Q,u) of N = 20, 5=1 VBS model. It is made of 
delta-function peaks at each excited state. The strength is expressed as percentage of total strength 
at the momentum Q. The delta-function peak enters into the continuum and changes to a round 
peak. Both are very close to the variational energy given by eqn.(12). 
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FIG. 9. Fraction of the total scattering intensity concentrated in the lowest energy state of 
S = 1/2 Heisenberg chain for N = 16,20,24 chains. The size-dependence is strong in all areas of 
momenta. This means that the lowest energy states are the bottom of the continuum made by two 
spinon excitations. 

FIG. 10. Energy of the first and second excitations of S = 1/2 isotropic Heisenberg point. It 
seems that there is no gap between them. 

FIG. 11. Structure factor S(Q) of S = 1/2 isotropic Heisenberg chain. It has logarithmically 
divergent peak at Q = tt and behaves as \Q\. 

FIG. 12. Dynamical structure factor S(Q,u>) of N = 24, S = 1/2 isotropic Heisenberg model. 
It is made of delta-function peaks at each excited state. The strength is expressed as a percentage 
of total strength at the momentum Q. 
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